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_ . Abstract 
(^ . 

We establish an indirect link between relic neutrinos and the dark energy sector which originates 

(~| , from the vacuum energy contributions of the neutrino quantum fields. Via renormalization group 

Qh' effects they induce a running of the cosmological constant with time which dynamically influences 

Q \ the evolution of the cosmic neutrino background. We demonstrate that the resulting reduction of 

5-H ' the relic neutrino abundance allows to largely evade current cosmological neutrino mass bounds 

c/3 \ and discuss how the scenario might be probed by the help of future large scale structure surveys 

, ^, ' and Planck data. 

(N 

^ ■ 1 Introduction 

O 



Two of the outstanding advances of the last decades were the confirmation of neutrino oscillations 
and the observation of a late-time acceleration of the cosmological expansion attributed to a form 
of dark energy, its simplest origin being a Cosmological Constant (CC). In this work, we discuss the 

C , implications for cosmological neutrino mass bounds arising from an interaction between the CC and 

^"^ ' relic neutrinos. 



X 



According to Big Bang theory the Universe was merely about one second old, when cosmic neutrinos 
fell out of equilibrium due to the freeze-out of the weak interactions. Prom the standard presumption 
C^ \ that the neutrinos have been non-interacting ever since decoupling, it follows that their distribution 
is frozen into a freely expanding Fermi-Dirac momentum distribution. Consequently, these neutrinos 
are assumed to homogenousely permeate the universe as the cosmic neutrino background (Cz/B) with 
a predicted average number density of ^ n^^^ = ^ np^ = 168 cm"'^ today. This substantial relic abun- 
dance only falls second to the photons of the cosmic microwave background (CMB), the exact analog 
of the Ci^B. However, owing to the feebleness of the weak interaction, all attempts to directly probe 
the Ci^B in a laboratory setting have so far been spoiled [H El [Sj I5- Yet, evidence for its existence 
and knowledge of its properties can indirectly be inferred from cosmological measurements sensitive 
to its presence like the abundance of light elements as well as CMB and Large Scale Structure (LSS) 
measurements (see e.g. Ref. p] for a recent review). For example, LSS data can reflect kinematical sig- 
natures of relic neutrinos which arise, since on scales below their free streaming scale massive neutrinos 
cannot contribute to the gravitational clustering of matter. Through the metric source term in the 
perturbed Einstein equations, this lack of neutrino clustering translates into a small scale suppression 



of the matter power spectrum which depends on the fractional energy density provided by neutrinos 
to the total matter density. Hence, within the standard cosmological model, LSS data allows to infer 
bounds on the absolute neutrino mass scale ^ rrii, = 0.2 -=- 1 eV {2a) depending on the employed 
data sets (see e.g. [U [7l El (U HOl dU [l2l [131 EH ES])- However, these bounds are in tension with the 
claim of part of the Heidelberg-Moscow collaboration of a > 4(T evidence for neutrinoless double beta 
decay translated into a total neutrino mass of ^ruy > 1.2 eV at 95% c.l. f\3\. While this claim is 
still considered as controversial (see e.g. [E]), at present other independent laboratory experiment 
are not sensitive enough to verify it. However, considering that cosmological neutrino mass bounds 
strongly rely on theoretical assumptions on the Ci^B properties such an apparent tension might hint 
at new physics beyond the standard model. Namely, it might imply the presence of non-standard 
neutrinos interactions which reduce the relic neutrino abundance and thus relax cosmological neutrino 
mass bounds (see Refs. [SI EH] and references therein). Hence, the improvement of current mass limit 
from tritium /3 decay experiments, J2''^v < 6.6 eV {2a) [El [20], by the approved tritium experiment 
KATRIN with a projected sensitivity of 0.2 eV [2TJ [22], offers the exciting possibility to probe such 
new neutrino interactions in the near future. What kind of new interactions could neutrinos have? 

Recently, it has been argued that a natural origin for new neutrino interactions would be the sector 
responsible for dark energy, a smooth energy component, which is assumed to drive the observed 
accelerated expansion of the universe [231 [Ml [25l [26] . A scenario has been proposed [21], in which the 
homogenousely distributed relic neutrinos are promoted to a natural dark energy candidate, if they 
interact through a new scalar force. As a clear and testable signature of this so-called Mass Varying 
Neutrino (MaVaN) Scenario the new neutrino interaction gives rise to a variation of neutrino masses 
with time. The resulting phenomenological consequences have been explored by many authors, in 
particular the implications for cosmology and astrophysics have been considered in Refs. [271 EH [29l 
[301 [an [31 [331 [Ml [35]. 

In this paper we propose an alternative approach for the realization of a dynamical influence of the dark 
sector on relic neutrinos (and vice versa), which is complementary to the avenue in the MaVaN Scenario 
in the following sense. While we also start from the requirement of energy-momentum conservation 
of the coupled two-component system, we, however, require the neutrino masses to be constant. As 
a direct consequence, the energy exchange with the dark sector demands a non-conservation of the 
number of neutrinos and accordingly its variation with time. Such a set-up naturally arises on the 
basis of Einstein's theory of general relativity in the appropriate framework provided by quantum field 
theory on curved space time [36] ■ Namely, this description requires the running of the vacuum (or zero 
point) energy or equivalently the CC as governed by the renormalization group [37l [38ll39l [40l [411 142] . 
which implies an energy exchange between the CC and matter fields [Ml HH HH [l6] . Within a suitable 
renormalization scheme, we find the relic neutrinos to govern the running of the CC at late times, when 
they turn non-relativistic. We analyse the consequences of the energy-momentum exchange between the 
relic neutrinos and the CC both for the late-time dynamics of the running CC and for the phase-space 
distribution and the abundance of relic neutrinos. Interestingly, independent of the absolute neutrino 
mass scale realized in nature, in the redshift range accessible by current cosmological probes we find 
the effective equation of state of the running CC to be extremely close to —1 as suggested by recent 
data (e.g. [3^). However, we demonstrate that within the considered scenario the effects of the Cz^B 
on cosmological measurements are drastically reduced. As a consequence, current neutrino mass limits 
arising from cosmology can be largely evaded. Furthermore, we discuss how the considered scenario 
can be tested in the near future by the help of cosmological and laboratory based measurements. 

The outline of the paper is as follows. In Sec. [2] we provide the theoretical framework for the considered 
scenario, and in Sec. [3] we analyse in detail the time evolution of the neutrinos and the CC arising from 
their interaction. Sec.[4|deals with the phenomenological consequences of the new neutrino interaction, 
in particular, with the implications for the inference of neutrino mass bounds from current LSS and 
CMB data. In Sec. [5] we provide an outlook how the considered scenario can be tested both by future 



cosmological probes and by experiments sensitive to Pauli-Blocking effects arising from the presence 
of the Cz^B. 



2 Running vacuum energy 

The existence of zero-point energies is an inevitable consequence of quantum field theory in which every 
quantized field is contributing to the vacuum energy density A. In the absence of notable gravitational 
interactions, however, the CC can be completely ignored as it has no influence on non-gravitational 
physics. In contrast to this, the cosmological evolution can be very sensitive to the CC |l8j, and 
the impact of quantum corrections should be taken into account. Unfortunately, the calculation of 
the vacuum contribution A A of a quantum field of mass m leads to a divergent result, which can be 
expressed as a momentum integral over all zero modes: AA ~ J d^p a/p^ + m? . Naive estimations 
with a high-energy cut-off, e.g. at the Planck scale Mpi, are far from successful since they lead to huge 
discrepancies with the measured value Aq of the CC: 

AA ~ M^i f» 10^2^ Aq. 



Nevertheless, quantum field theory on a curved spacetime [36] allows us to obtain at least some infor- 
mation from the divergent contributions. In this framework gravity is considered to be a completely 
classical theory well below the Planck scale, which provides the background for quantum fields and 
their excitations. By using well-known procedures, infinities like AA arising in this set-up can be made 
finite and eventually lead to meaningful results by treating them in an analogous way as for example 
divergent quantum corrections of the fine-structure constant in quantum electrodynamics (QED). Con- 
sequently, regularization and renormalization techniques allow to formulate a renormalization group 
equation (RGE) for the CC. Of course, this procedure does not determine the absolute value of A, but 
only its change with respect to an initial value Aq that has to be measured by experiments. Similar 
to the running fine-structure constant, the vacuum energy A(/x) becomes a scale-dependent variable 
with /i as a characteristic energy scale, called the renormalization scaled For each bosonic/fermionic 
field degree of freedom (DOF) the RGE can be expressed in the form 

C^A Tft 

where m denotes the mass of the corresponding quantum field. The complete equation is then given by 
the sum of DOFs of all fields in the theory. For free fields Eq. ([T| is an exact result, and even in curved 
space-time there are no higher order corrections. Unlike QED, where the renormalization scale fi is 
usually related to external momenta or energy scales, the running of A(/i) originates from a zero-point 
function and does not involve external scales. Therefore, we have to find a suitable identification for /i. 
In cosmology several choices have been discussed before in the literature, e.g. the Hubble expansion 
rate H or quantities related to horizons, see [HI EOl EI] and the references therein. However, in 
contrast to such global choices for /i, in this work we investigate a scale that is more characteristic for 
the contributing field and the corresponding particle spectrum. To motive its choice, it is important 
to note that Eq. ([T]) has been derived in the MS renormalization scheme that usually does not show 
decoupling effects at energies sufficiently below the mass m [52l [53l [54] • Since the current CC is 
much smaller than m^ for most particle masses, the naive evolution of the RG running down to low 
energies fi <^ m might induce unrealistically huge changes in the CC at late times. Instead of stopping 
the RG evolution by hand in the low energy regime, we choose our RG scale such that it naturally 
implements some sort of decoupling when fj, approaches the field mass m. 



^One also finds a RGE for Newton's constant G, but its running is suppressed by the ratio (m/AIpi)'^. Therefore we 
can safely ignore this scale-dependence in this work. 



Before identifying the renormalization scale, we have to note first that a time-dependent // generally 
implies that the vacuum energy A cannot remain constant anymore, even though its fundamental 
equation of state is still —1. Therefore, in order to have a consistent framework, Bianchi's identity has 
to be satisfied. On a cosmological background this means 

A + p + 3H{p + P) = 0, (2) 

where the dot denotes the derivative with respect to cosmic time t. For concreteness, let p and P be 
the total energy density and pressure of all components apart from the CC. Since for self- conserved 
matter the term p + 3H{p + P) vanishes, there must obviously exist a non-trivial interaction between A 
and p accounting for the exchange of energy-momentum between both sectors. Consequently, both the 
matter energy density and particle number density do not obey the standard dilution laws anymore. 
Note that this is qualitatively different from self-conserved dark energy sources, like (non-interacting) 
scalar field or holographic dark energy models. 

Now, one might ask which part of the total matter energy density consisting mainly of cold dark 
matter, baryons, leptons and photons is actually exchanging energy-momentum with the running CC 
to compensate for the time dependence of A in Eq. ([2]). Since the microscopic realization of the implied 
interaction is not known, let us in the following motivate our approach to answer this question. The 
vacuum contribution of a quantum field of mass m leads to a RGE as specified in Eq. ^. At the same 
time, the quantum field describes particles of mass m, whose energy density contributes to the total 
matter density p. Accordingly, a natural way of satisfying Eq. Q seems to be provided if the system of 
the vacuum contribution and the energy-momentum of the quanta arising from the same quantum field 
is conserved. Correspondingly, each field obeys an equation of the form of Eq. ^, however, with A, 
p and P belonging to this field only. Following the same line of reasoning, the RG scale p should be 
related to this self- conserved system, too. 

To examine the RG evolution of a weakly interacting quantum field, in our framework it is appropriate 
to introduce the momentum distribution function /, which characterizes the density of the correspond- 
ing particles in a given momentum bin on the homogeneous background of the universe. By its help, 
the energy density, pressure and the particle number density of each field DOF can respectively be 
expressed as 



n 



——^jd?yu){y)f{y), uj{y) := '^ y"^ + j^ , 

Here, we take / to depend on y := p/T, the ratio of the momentum p and the effective temperature T 
characterizing the species, which in the considered framework have the same scaling with redshift such 
that y is constant in time. Hence, while the effective temperature T is cooled by the expansion, the 
form of such a distribution is maintained in the absence of interactions. However, as it will turn out 
later on, taking into account the interaction with the CC, the distribution function will not expand 
freely anymore but also change with time. 

Accordingly, a reasonable choice of the RG scale seems to be the average energy p/n of the system, since 
it behaves at high energies like the characteristic momentum or temperature scale T. Furthermore, it 
is bounded from below by the field mass ttt., thereby implying the end of the RG running. In order to 
keep the spin statistics intact, we also implement the factors (1 ± /) for bosonic/fermionic fields. 



fi (27r)^ 



j d?yujf{l±f), n:=-^ j d?yf{l±f). (4) 






Note that by these means p represents the average energy-density of the system which is available for 
the interaction with the CC. Furthermore, we would like to point out that for more general systems or 
sets of individual particles, the resulting form of Bianchi's identity and the RG scale might look very 
different from the considered case. 



3 The RG evolution of A and the Cz/B 

In the following, we restrict our discussion to homogeneously distributed matter in the universe, that 
is to CMB photons, to background neutrinos and in principle to possibly existing other relics. The 
photons in the CMB are massless and thus not affected by RG effects. However, relic neutrinos, 
should be subject to RG running since they are massive and very weakly interacting. According to Big 
Bang theory, their distribution was frozen into a freely expanding, ultra-relativistic Fermi-Dirac (FD) 
form, /o(y) = l/(exp(y) + 1), when neutrinos decoupled from the thermal bath at a temperature scale 
of Tdec ~ 1 MeV ^ m. In the following, we will therefore use this form as initial condition for / and 
calculate its time evolution by virtue of Bianchi's identity ([2|) and the integrated form of the RGE ^ 

AM = A„-^l„^^ (5) 

Let us stress that the following considerations generically hold for any weakly interacting fermionic 
species of mass ra exchanging energy-momentum with the CC, if its initial distribution function assumes 
the ultra-relativistic FD form fo{y)- For this reason, and for simplicity, we will omit the index u in the 
following, but reinstate it when we investigate the phenomenology of interacting neutrinos in Sec. HI 
Accordingly, without loss of generality, we have taken both p, n and P in Eqs. ([3|) as well as A(//) to 
correspond to one fermionic DOF, since the actual number of fermionic DOF drops out in Eq. ([2]). 
Here, the RG running of the mass can be neglected due to the absence of substantial interactions with 
other fields. 

Because of the momentum integrals in Eqs. (I3|), the Bianchi identity ([2]) yields a rather complicated 
non-linear integral equation for the time-evolution of f{y): 



= p + A + 3H{p + P) = ^L-Jd^yco 
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327r^ '' ^ "" ' p \ Tujfi 

While the interaction tends to deform / away from its equilibrium form /o as discussed in the following, 
we will still use the effective temperature T = To(2: -|- 1) as a measure for cosmic redshift as in 
the standard case. Accordingly, the present characteristic temperature of FD distributed neutrinos 
is To ^ 1.7 ■ 10~^ eV. To find a solution of the integral equation, we require the expression under the 
integral to vanish for any momentum value y. This removes the outermost integral, and we obtain an 
equation that can be solved numerically. In terms of the derivative /' with respect to the redshifcl z 



V and /' are related by f{y) = -{z + l)Hf'{y). 



we find the following evolution equation 



8 A^ 

,2 A 



N ■■= ^ + 1-7^, (7) 



_r 

^ := — ^ + (1 - 2/) f 7^ - iV (8) 

Let us in the following analyse this equation and compare it with our numerical results for the evolution 
of /. First, as a consequence of employing p instead of p in Eq. (|1|, one can see that the factor /(I — /) 
on the right-hand side ensures Pauli's principle to hold, < / < 1. Moreover, at early times, when p ~ 
r^ » 771^, the dynamics is strongly suppressed, |/'[ ^ 1, and the neutrinos just behave as in the 
uncoupled case. During this time the denominator D is positive for any value of y, while both the sign 
of the nominator A'^ and of /' are y-dependent. For large momenta, y ^ oo, we find f',N > implying 
the decay of neutrinos. Considering that p, is decreasing, this is the expected behavior following from 
Eq. (121). However, in the small momentum regime, y — > 0, the distribution function increases because 
of /', A^ < 0. Since we are lacking a microscopic description of the dynamics of the coupled system, 
just from solving Bianchi's identity it is impossible to say whether the growing / owes to slowing down 
high-y neutrinos or is also due to particle production. 

Finally, the distribution function is either approaching / = or / = 1, when the positive denominator in 
Eq. do]) becomes zero. More precisely, the first zero appears for large momentum, where /(y —i- oo) — > 
and w — > (X). Since the second term in Eq. ([8]) in this limit is —1, the time when this happens is 

4 

determined by p = ^^- This roughly corresponds to T k, m/3, i.e. the time when neutrinos turn 
non-relativistic. To arrive at this estimate, we have assumed that / has not changed much yet from its 
initial FD form /o and in addition that /> ~ p. At the momentum position of this zero the spectrum is 
simply cut off with all larger momenta being removed, which is illustrated nicely in the last three plots 
in Fig. [U for decreasing redshift. This figure also demonstrates that for increasing time and further 
decreasing /), the momentum cut-off is moving to lower y- values. 

A second zero in the denominator D will eventually occur for small values of y. This happens when p 
has decreased sufficiently and (1 — 2/) is negative. In the corresponding y-range we find f'.,N < 0, 
and accordingly in this case / acquires its maximal value 1. As illustrated in the last of the two plots 
in Fig. [H with decreasing redshift, the second zero of D occurs for larger and larger y- values for which 
consequently /(y) — > 1. Finally, both zeros in D meet at y = i/ei ~ 0.7, which implies that / cannot 
be evolved anymore according to Eq. ©. 

At this stage the spectrum represents a maximally degenerate Fermi-Dirac distribution with all states 
completely occupied up to the Fermi momentum i/ei- Once this Fermi sphere has formed, it is still 
possible to pursue the further time evolution of /(y) analytically by allowing the Fermi momentum ye(t) 
to be time-dependent. As a result, we find respectively for the RG scale and the neutrino energy density 

P = ^^VymT' + m\ p = —^ / dyyV (9) 

n 27r^ Jq 

Accordingly, other /(y)-dependent quantities can be determined by replacing the distribution function 
by a step function which is 1 for y < ye{t) and elsewhere. From Bianchi's identity we therefore obtain 
a much simpler differential equation for ye{t), 

_ Hy^jt) f >>. Hy^jt) 

'^^^"l-16ye(t)S..3 - l-16X,,(t)- ^''^ 
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Figure 1: Evolution of the momentum distribution profile f{y) of one neutrino with two DOFs and 
mass nil, = 9 x 10~^ eV. Note the different scaling of the axis of ordinates and of the abscissae at each 
redshift z. The plot corresponding to z = 7.12 demonstrates that the denominator D{y) in Eq. ([8]) 
exhibits two zeros, at which f{y) becomes 1 or 0, respectively. The final profile is acquired at z r^ 4.75. 
As indicated by the shaded region in the last z = 5.47 plot, it corresponds to a maximally degenerate 
Fermi sphere with all states occupied up to the Fermi momentum ye ~ 0.57 (dashed-dotted line). 



As can be read-off from Eq. (fTOl) . since neutrinos are already non-relativistic and thus y ^ Hi Ueit) can 
only be evolved until it takes the final value yeit) ~ m/(16T). In our system this means that neutrinos 
are not able to exchange energy-momentum with the CC anymore. This is a direct consequence 
of including the Pauli principle in the RG scale, which blocks any further change in the coupled 
system. Therefore, it can be understood as the decoupling of the neutrinos from the running CC. 
Thus, henceforth, both the running of /i and of the CC is taken to cease such that the neutrinos evolve 
independently. 

To find a (semi-)analytical expression for the final Fermi momentum ye, in the following, we will 
estimate the redshift Tg = To(ze -|- 1) corresponding to the decoupling time. To this end, we proceed 
by approximating the redshift evolution of interacting neutrinos. Accordingly, when the interaction 
effectively sets in, the temperature roughly satisfies Ta = m/3, and we approximate the distribution 
function f{y) by a Fermi sphere with initial Fermi momentum ya > ye- At this time, ya can be 
determined by comparing the energy density psph of this sphere with the energy density ppD of a 
standard FD distribution. In the non-relativistic neutrino regime, we obtain from the equality of 



1 ^3 s , 3C(3) o 

Psph^^T^my^ and ppD ^ ^;^T^m, 

the initial Fermi momentum ya which takes the form ya = (|C(3))^'^ 



(11) 
1.76. By integrating the 



differential equation p^ in the formal 



2; + l l-16-ye(z) 

from Ta = m,/3 to Tg with y^ = m/(16Te), we find ye = yaQxpiTn/{16yaTa) — 1) ^ 0.72 and thus Tg « 
m/11.5. Accordingly, the estimate for Tg is of the same order of magnitude as the numerical result Tg « 
m/9. Note that the final value, 

ye = m/{l6Te) ^ ^, (13) 

turns out to be independent of the mass parameter, which is confirmed by our numerical calculation. 

Because the maximal physical momentum pe = yeTe = -j^ of the neutrino spectrum is much smaller 
than the mass m, after decoupling from the CC, the neutrinos redshift just like non-relativistic dust, 
thus contributing to the dark matter. However, since they have transferred a substantial amount of 
energy-momentum to the vacuum energy during the interacting phase, with respect to non-interacting 
neutrinos of equal mass and energy density /OpD, the magnitude of their energy density is reduced by 
a factor, 

^^ = = ° flr\m ^ 30 for T < Te. 14 

In the considered scenario, this relation also defines the apparent neutrino mass rriapp, which follows 
from assuming the standard scaling law for the number density npo and equal energy density, p ~ 

W-app'^FD for T <Te. 

Let us in the following qualitatively discuss the impact of the coupling between neutrinos and the CC 
on the background evolution. From the previous discussion we first note that the main effects of the 
interaction occur during an epoch set by the neutrino mass, 

Te^-<T<Ta^-. (15) 

Accordingly, the smaller the neutrino mass the later does the energy-momentum exchange between 
the coupled system occur. However, also the impact on the background evolution decreases with 
the neutrino mass, while the running of A becomes entirely negligible for m <C Aq according to 
Eq. ([1|). In addition, in comparison to the ACDM model with massless neutrinos, the energy density 
in light massive neutrinos is not much larger. As an example, we have plotted the redshift evolution 
of the cosmological constant in Fig. [2] for a hierarchical spectrum with one massless and two massive 
neutrinos as allowed by neutrino oscillation experiments. At high z both of the massive neutrinos 
contribute to the running of A, while the heaviest neutrino of mass m^ = 5 x 10~^ eV is the dominant 
one. At z ~ 30.6 its momentum distribution has acquired its final, maximally degenerate form and 
thus the neutrino decouples. Afterwards, the second lightest neutrino with mass m2 = 9 x 10~^ eV 
is left over to drive the CC to its current value Aq, which is reached at z Ri 4.75. Neutrinos lighter 
than To/9 th 1.5 x 10^'^ eV could be interacting with the CC even today, but their influence on the 
running A would be too small to be detected. 

However, we encounter a more pronounced effect of the interaction for larger masses, where A was very 
negative in the past and evolves to its current observed value Aq > 0. A negative CC might first look 
strange, but there are many potential sources for negative vacuum energy, e.g. spontaneous symmetry 



3 



The general solution of Eq. I|12p for arbitrary j/a.e and Ta,^ reads 

fmfl 1 

ye = ya exp — — 

\16 yiaVa J-eVe 



SI 




1+z 

Figure 2: Vacuum energy density A.{z)/Aq as a function of redshift for a hierarchical neutrino mass 
spectrum with mi = eV, m2 = 9 x 10~^ eV, ma = 5 x 10~^ eV and 2 DOFs each. Note the different 
scahng of the ordinate. 



breaking [39] or supersymmetry [55j. In addition, it turns out to be harmless since in the relativistic 
regime, where T ^ m and /u ~ T the vacuum contribution is given by 
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Obviously, this is much smaller than the energy density ~ T^ of the dominant radiation components 
at early times. Therefore, the decay of the neutrinos and the influence of the running CC on the 
background evolution is negligible during this stage, thus leaving Big Bang Nucleosynthesis (BBN) 
unaffected by this mechanism. Only in the redshift range, where the energy exchange with the neutrinos 
is relevant, A may have a larger impact on the background evolution. To illustrate the effect we consider 
the case of three degenerated neutrinos with mass rrii, = 2.2 eV and 2 DOFs each, as maximally allowed 
by /3-decay experiments. 

In the redshift range z ~ 1400 -r- 3900 a large fraction of neutrino energy-momentum gets transferred 
to the CC and thus the acceleration factor q = ^ indicates a prolonging of the radiation dominated 
epoch (g « — 1) with respect to standard cosmology. In other words, temporarily, the expansion is 
decelerating to a greater extend, as shown in Fig. [31;). The reason for this background behavior is that 
the coupled negative vacuum energy and the decaying neutrinos act as a very stiff system as illustrated 
by the corresponding effective equation of state Weg in Fig. [3^). In the final stages of the interaction 
the universe quickly approaches the standard evolution of a A Mixed DM model (AMDM) again, which 
corresponds to a ACDM model including neutrino masses. Therefore, the transition from the radiation 
dominated to the matter dominated epoch is much sharper than in the non-interacting case. As of 
then, the CC exhibits its standard value and redshifts with an equation of state of a; = —1, while the 
neutrino energy density dilutes according to the standard law, however, its magnitude is reduced by a 
factor 30. Note that in the rrii, = 2.2 eV case this happens before recombination. 

Our third example shows neutrinos with degenerate masses rriy = 0.5 eV. In this case, the energy- 
momentum transfer starts after the time of recombination at approximately z = 900, and in comparison 
to the case of rriu = 2.2 eV has less impact on the background evolution due to the smaller neutrino 
mass. Nevertheless, Fig. [3] suggests that the integrated effects on the expansion rate from today to 
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Figure 3: Cosmological evolution in the case of three degenerate neutrinos with two DOFs and 
mass mu = 2.2 eV (left) and rrii, = 0.5 eV (right) each, a) The relative energy densities Qi of the com- 
bined neutrino-CC system (solid), neutrinos (short-dashed), the CC (dashed-dotted), matter (long- 
dashed), radiation (dotted), respectively, b) The effective EOS Wgff of the combined neutrinos-CC 
system (solid) in the interacting phase, the EOS of neutrinos (dashed) and the CC (dashed-dotted), 
respectively, c) The acceleration factor g = f| within ACDM (dotted), AMDM (dashed) and AMDM 
with neutrino-CC interaction (solid), respectively, d) The ratio of the Hubble expansion rate within 
AMDM with neutrino-CC interaction and ACDM (solid) as well as AMDM and ACDM (dotted), 
respectively. 

the time of recombination have an in principle measurable influence on cosmological observations, cf. 
SecH 

The results discussed so far should also be applicable for other fields of higher masses as long as they 
are not interacting too strongly. Then, the mechanism starts at earlier times, when A was even more 
negative. However, as argued above, we expect the CC to be sub-dominant also during these stages 
and thus not to seriously interfere with standard Big Bang cosmology. Let us justify this statement 
by explicitly calculating the maximum of the ratio between the shift AA of the running CC and the 
energy density p of the interacting fermions. 



AA 



r := 



m 



32 \ye{t)T 



In 



ye{t)T 



fJ-e 



(16) 



Note that the number of fermionic DOFs drops out. Since the maximum occurs close to the end of the 
interacting phase, we have employed the energy density and RG scale of a Fermi sphere from Eqs. (fTT]l 



and ([9|). In addition, the final scale p^ = |||"i^ is fixed by VeTe = m/16. Therefore, the maximum 
of r is acquired for yT/m ^ 0.11, and its value r^ax ~ 0.57 turns out to be independent of the field 
mass. Thus, let us apply this upper bound to the early universe. Accordingly, deep in the radiation 
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dominated epoch, the total energy density ptot is dominated by many relativistic particle species such 
as photons, neutrinos, electrons and so on. In case one of these species becomes non-relativistic, it 
induces a vacuum shift which maximally corresponds to a small fraction of ptot- Later on, however, 
when only few relativistic species are left over, this fraction can be as large as r^a.x- Therefore, the 
running vacuum energy shows a tracking behavior, which becomes better and better as time advances. 

Despite the smallness of r, one might still worry about early dark energy constraints at BBN. At this 
time, electrons and positrons should be subject to our mechanism and the induced changes in the CC 
could harm the background evolution. However, from particle species like e^ that annihilate after 
their decoupling from the thermal bath, there is not much energy-momentum left over that can be 
transferred to the CC. Therefore, the interaction of the particles with the CC effectively ceases at this 
moment, and the induced vacuum shift is much smaller compared to a fully developed mechanism. 

It is remarkable that even though we cannot explain the smallness of the current CC, we have found 
a mechanism that keeps its quantum corrections under control during the cosmological evolution. 
Additionally, in contrast to models of tracking quintessence, the considered scenario only requires a 
single free parameter, Aq. It is also worth noting that the final vacuum shift induced by the neutrinos 
is of comparable magnitude as the present CC value, although both quantities are apparently not 
directly related to each other. 



4 Phenomenological implications 

4.1 Relaxing cosmological neutrino mass bounds 

In this section, we will consider the consequences arising from the non-standard neutrino interaction 
with the CC for the inference of neutrino mass bounds from current as well as from future cosmological 
measurements such as of the CMB and of LSS. 

We will start by extending the discussion of the last section on modifications to the background 
evolution caused by the interaction between the neutrinos and the CC. This will allow us to give 
a qualitative estimate of its impact on the spectrum of the CMB fluctuations and of the arising 
possible restrictions on the masses of interacting neutrinos. Furthermore, it will set the stage for 
our subsequent (semi-) analytical analysis in Sec. 14. 1.2) which demonstrates that within the considered 
scenario all current bounds on neutrino masses derived from LSS measurements can be evaded. A 
detailed numerical analysis can be found elsewhere |56| . 



4.1.1 Relaxing neutrino mass bounds from CMB data 

As already mentioned in the last section, the background evolution is modified in the presence of the 
considered non-standard neutrinos interaction due to the following two colluding effects. The first 
one results from the stiffness of the effective equation of state of the neutrino - CC system during the 
interacting phase which tends to prolong the radiation dominated regime as argued before. The second 
effect on the expansion rate arises generally if the ACDM model is extended to include neutrino masses 
(AMDM) [5]. However, as discussed in the following, its impact turns out to be much less important 
in the presence of the considered non-standard interaction. 

In general, while massive neutrinos are still relativistic at sufficiently early times and are thus counted 
as radiation, after the non-relativistic transition they contribute to the dark matter. This implies that 
the relative density of other species or the spatial curvature today have to get modified in comparison to 
a universe with massless neutrinos. However, for interacting neutrino the energy-momentum transfer 
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to the CC becomes efficient as soon as they turn non-relativistic. This ensures that within a relatively 
short phase of interaction the neutrino energy density is reduced to 1/30 of its standard value, cf. 
Eq. ([HI). Accordingly, in contrast to the standard case, the contribution of interacting neutrinos to 
the total energy density today is Qi^h'^ < 2.4 x 10~^ (for ^ tji^, < 6.6 eV as required by tritium 
bounds [l9j[20]). Hence, for a flat universe, this demands a negligible reduction of the dark matter 
density compared to the massless neutrino case. As a consequence, the usual effect of neutrino masses 
on cosmology of postponing the matter-radiation equality is much less pronounced for interacting 
neutrinos. However, the interaction with the CC has an influence on the background evolution and 
thus on the time of transition from radiation to matter domination, where its impact grows with 
increasing neutrino mass. This can be seen by comparing the evolution of the acceleration parameter 
q with and without the interaction in Fig. [3l While the radiation dominated regime is prolonged, the 
transition to the matter-dominated regime happens faster in the presence of the interaction, since it 
quickly causes neutrinos to become highly non-relativistic. This is also demonstrated by Fig. [Sj where 
the ratio of the Hubble expansion rate within AMDM with neutrino-CC interaction and ACDM with 
massless neutrinos is compared to the corresponding ratio within AMDM without the new interaction 
and ACDM. Especially the effects on the integrated expansion rate are expected to be smaller for 
interacting neutrinos as suggested by Fig. [3l For the CMB power spectrum, however, it is mainly the 
integrated expansion rate, which determines the important scales that set the location and the height 
of the acoustic peaks [5]. More precisely, the sound horizon at recombination is given by, 

oo 

■Srcc 

with Cs denoting the sound speed of the baryon-photon fluid. The observable angular scale of the 
acoustic oscillations depends on an integral from the present to recombination involving the density 
and the equation of state of the running CC [57l[58]. While keeping all other parameters fixed, with 
respect to a ACDM model, in principle a linear shift in the angular position of the CMB features could 
be induced. However, the effective ujes only differs from u; = — 1 during the relatively short phase of 
interaction to an extent which decreases with the neutrino mass. Furthermore, already much before the 
dark energy dominated regime the energy density in the running CC has essentially reached its present 
value (cf. Fig. [2]). Accordingly, the integrated result generally will not differ much from the ACDM 
caseO We note that secondary effects on the CMB acoustic peaks resulting from the free streaming 
of relativistic neutrinos (i.e. a phase shift and an amplitude reduction at large multipoles) [59\ are 
expected to essentially be unaltered in the considered scenario, since in the relativistic neutrino regime 
the coupling between neutrinos and the CC is very strongly suppressed (cf. Sec.[2l). 

Hence, in summary, qualitatively, compared to the standard ACDM case with massless neutrinos, the 
modifications to the CMB spectrum arising from the inclusion of neutrino masses are less pronounced in 
the presence of the non-standard neutrino interaction. Therefore, for interacting neutrinos the current 
neutrino mass bound of ^ rrii, = 2 -^ 3 eV (at 95% c.l.) gained from WMAP3 data alone is expected to 
be relaxed (cf. [5] and references therein, see [60] for the tightest bound). Thus, a priori, the considered 
scenario does not seem to be excluded by current CMB data, even if the sum of neutrino masses is 
in the super-eV range as allowed by current tritium experiments. A detailed analysis of the resulting 
modifications to the CMB spectrum and the arising bounds on the masses of interacting neutrinos goes 
beyond the scope of this paper and will be investigated elsewhere [56]. 

In the next subsection we will discuss the consequences of the non-standard neutrino interaction on 
the inference of neutrino mass bounds from LSS measurements. 



■^Prorn similar arguments follows that the influence of the interaction on the magnitude of the integrated Sachs- Wolfe 
effect at low multipoles is expected to be small. 
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4.1.2 Evading current neutrino mass bounds from LSS measurements 

Structure formation is sensitive to neutrino masses through kinematic effects caused by the neutrino 
free streaming as characteristic for hot and warm dark matter [6H Il2]. The associated important 
length scale corresponds to the typical distance neutrinos can propagate in time t in the background 
spacetime and is governed by their mean velocity (v) , 

where Aps/o and fcps respectively denote the comoving free streaming length and wavenumber of the 
neutrinos. As a consequence, on the level of perturbations neutrinos do not aid to the gravitational 
clustering of matter on scales below the horizon when they turn non-relativistic. However, their energy 
density py contributes to the homogenous background expansion through the Priedmann equation. 
Through the metric source term in the perturbed Einstein equation, this imbalance leads to a slow 
down of the growth of matter perturbations on small scales (cf. Ref. [5] for a recent review). It crucially 
depends on the relative fraction provided by neutrinos to the total energy density p^ of mattei|j, 

Pm 15 eV 
where pm comprises the energy densities of cold dark matter, baryons and neutrinos. 

Assuming f,/ ^ 1, on scales smaller than the horizon when neutrinos turn non-relativistic, the resulting 
net effect on the present day matter power spectrum Pmik) for a normalization at k — > is a 
suppression of [62] 

In the following we will demonstrate that the proposed interaction between relic neutrinos and the 
CC allows to completely evade present neutrino mass limits deduced from LSS data below the tritium 
bound nil, < 2.2 eV [HJ [20], because the neutrino energy density gets reduced below the current 
sensitivity of LSS measurements on all accessible scales. However, at the end of the section we will 
argue that both for the Planck mission as well as for weak lensing or high redshift galaxy surveys it 
seems feasible to probe the proposed non-standard neutrino interaction. To this end, we will investigate 
the modifications to the characteristic scale for free streaming signatures induced by the considered 
neutrino interaction. 

Let us start by noting that for structure formation the interaction generally has a negligible influence 
on the expansion rate independent of the actual neutrino mass scale realized in nature. Namely, if the 
interaction takes place at late times during structure formation, the corresponding neutrino masses 
and the amount of energy-momentum transferred to the CC are too small to lead to an observable 
effect on H. However, in case the interaction has an impact on the expansion rate for sufficiently large 
neutrino masses, then it ceases already much before structure formation. Accordingly, in the following 
analysis we will adopt H{z) ~ H{z), where a tilde throughout this work labels quantities which assume 
an interaction between neutrinos and the CC. 

Let us now proceed by comparing the fractions fi, and f,y provided by neutrinos to the total matter 
density with and without non-standard interaction in the non-relativistic neutrino regime, T^, < m,^/3, 



'Note that the formula in Eq. lfT9|l approximately also holds true in case a very light neutrino state has not turned 
non-relativistic yet, because in this case its relative contribution to the energy density is negligible anyway [5J. 
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where Pu{z) ~ 7n^n,^{z). In accordance with Eq. ([HI, we define the apparent mass of the neutrinos 
also for times before the interaction ends, 

m^PpiTu) = m^^—— = m^ — — — . 21 

Pu[T,y) ny{Ti,) 

Here, rriapp represents the mass of a non-relativistic neutrino an observer would infer from the extent 
of a neutrino-induced suppression of the matter power spectrum, if he assumes the standard scaling 
of py for a Fermi-Dirac distribution. Per definition TTiapp absorbs the additional non-standard redshift 
dependence common to hy and py caused by the interaction with the CC. We have plotted mapp as a 
function of the redshift for different possible values of the neutrino mass in Fig. 3], which demonstrates 
that rriapp becomes constant for Ty < niy/^ when the energy exchange with the CC has ceased. Cor- 
respondingly, after neutrinos have effectively decoupled from the CC, their energy density again obeys 
the standard scaling law with redshift, py cc {1 + zY . However, compared to the energy density of a 
non-interacting neutrino of mass niy^ its absolute magnitude is reduced by a (mass independent) factor 
of TTi-iy/TTiapp — 1/30 (cf. Eq. (f2T]) and Fig. |4|. Consequently, since tritium /? decay experiments con- 
strain neutrino masses to be smaller than 2.2 eV [l9l[20], the apparent neutrino mass has to be smaller 
than TTT-app < 2.2/30 eV ~ 7 x 10"^ eV for Ty < my/9. Thus, according to Eq. (fT9l) . the suppression of 
the matter power spectrum on scales below the neutrino free streaming becomes proportional to, 

fy = ^^f. - ^fy - ^ < 0.015. (22) 

niy 30 450 eV 

This result demonstrates that for ^ ruy < 6.6 eV the interaction with the CC reduces the efficiency of 
the suppression of the power spectrum achieved by neutrino free streaming below the current sensitivity 
of LSS measurements. It should be noted that this result is independent of the actual scales where 
neutrino free streaming is relevant. It thus implies that within the considered scenario all current 
bounds on niy from LSS data can be evaded. 

However, in the next subsection, we will discuss the characteristic free streaming signatures of inter- 
acting neutrinos and see how they can possibly be revealed by future cosmological probes. 

4.2 Characteristic free streaming signatures 

Future weak lensing and also high redshift galaxy surveys in combination with Planck data promise a 
considerable increase in the sensitivity to neutrino mass to cri^my) < 0.05 eV [Ml El]. For instance, 
surveys of weak gravitational lensing of distant galaxies directly probe the matter distribution without 
having to rely on assumptions about the luminous versus dark matter bias in contrast to conventional 
galaxy redshift surveys. The improvement of the sensitivity is largely due to the possibility of break- 
ing parameter degeneracies in the matter power spectrum in particular by the help of tomographic 
information on the power spectrum resulting from a binning of the source galaxies by redshift. Accord- 
ingly, while existing LSS surveys are mainly sensitive to the transition region in the power spectrum, 
where free streaming effects start to become important, these future probes at different redshifts can 
accurately measure a broader interval of wavenumbers extending into the present non-linear regime 
(e.g. at the peak of the sensitivity of weak lensing surveys A: ~ 1 -^ lO/iMpc"^ at z = 0.5 [63], while 
galaxy surveys probe scales 4.5 x 10~'^/iMpc~^ -^ 1.5hMpc~^ at 0.5 < z < 6 [S^). By these means, 
they are in principle sensitive to the characteristic step-like effects on the power spectrum arising from 
neutrino free streaming and as a result to much smaller ^y. As we will see in the following, this offers 
the possibility to test the considered scenario by these future probes. The reason is that the interac- 
tion between neutrinos and the CC turns out to not only reduce the energy density, but also the free 
streaming scale of the neutrinos for a given mass. Accordingly, the free streaming signatures tend to 
be shifted towards the non-linear regime as we will demonstrate in the following. 
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To this end, let us proceed by comparing the characteristic velocities of neutrinos which determine 
the free streaming scale according to Eq. (fTSjl with and without non-standard neutrino interaction. 
Starting with the standard case, the neutrinos are assumed to be Fermi-Dirac distributed, the average 
neutrino velocity corresponds to the thermal velocity (v), 

iv)^^^ ^, (23) 

in,, rriy 

where {p) denotes the mean momentum of freely propagating relic neutrinos. However, as described in 
the last section, the presence of an interaction with the CC tends to remove high momenta from the 
neutrino spectrum. As a consequence, in this case the average neutrino velocity gets reduced and the 
transition to the non-relativistic regime is speeded up compared to non-interacting neutrinos. Since a 
detectable suppression of the power spectrum by future surveys corresponding to ^ m-app > 0.03 eV 
for the most ambitious projects [65l[5] requires the sum of neutrino masses in the interacting case still 
to be in the super-eV range, we can safely assume the interaction with the CC has ceased already much 
before structure formation. Thus, after the final stage of the interaction, the mean neutrino velocity 
{v) of the maximally degenerated neutrino distribution with Fermi-momentum p-p = tjeTu = 4/3 (p) 
can be expressed as, 

(^) - ^ - 7^. (24) 

Accordingly, by the help of Eq. (fTSl) and Eq. (fT3]l we arrive at the relation. 



iv) 64 
k-FS = ^FSTTT =i — A;fs - 7.11 /cFS, (25) 

(w) 9 

where H{z) ~ H{z) was assumed as justified above. Consequently, since on average the interacting 
neutrinos are slower, their comoving free streaming wavenumber k-p^ according to Eq. (|25p is roughly 
an order of magnitude larger than in the standard case. We would like to stress that this result holds 
independently of the actual neutrino mass scale realized in nature for redshifts 1 + z < 171,^/(9X0). 

For neutrinos becoming non-relativistic during matter domination, the comoving free streaming wave- 
number passes through a minimum at the time of the non-relativistic transition. Thus, in the standard 
case, the matter power spectrum is suppressed on all scales much larger than [5], 

/fcnr ^ 1.8 X 10-2 nll^ (^) '^' hMpC-\ (26) 

In contrast, the free streaming signatures of interacting neutrinos are expected to influence all wavenum- 
bers larger than the minimal comoving wavenumber at T^ ~ m,^/9, 

k^, ^0.23nl/' [^y hMpc-\ (27) 

For comparison, for Qm ^ 0.3 and 0.046 < m^y < 2.2 eV this corresponds to 2.1 x lO^^hMpc^^ < 
knr < 1.5 X 10-2/iMpc"^ in contrast to 2.7 x lO-^/iMpc"^ < Kr < 0.18/iMpc"^ Clearly, the free 
streaming signatures of interacting neutrino are shifted towards the non-linear regime. 

Hence, in summary, a clear signature for the considered scenario is provided by the non-standard 
correspondence between the extent of the suppression of the matter power spectrum due to neutrino 
free streaming and the corresponding scales which are affected. In other words interacting neutrinos 
appear lighter, however, at the same time can free stream a shorter distance in a Hubble time than in 
the standard case. 

Before turning to the prospects of testing the considered scenario by future probes, let us mention 
another characteristic signature expected to arise in its framework if ^ rrii, is in the super-eV range. 
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As argued above, in this case the new neutrino interaction has an impact on the time of matter-radiation 
transition whose extent grows with increasing neutrino mass. Accordingly, this leads to a translation 
of the turning point in the matter power spectrum in comparison to the ACDM case, because the 
scales entering the horizon at equality have a different size. However, the evolution of the acceleration 
parameter q is characteristically altered compared to the case of non-interacting massive neutrino as 
described above. Thus, it is expected to lead to a shift of the turning point compared to both the 
ACDM and the AMDM case whose location, however, has to be determined numerically [5^ . 

5 Probing the scenario 

5.1 Future cosmological probes and interacting neutrinos 

In this subsection, we will give an outlook how combined future weak lensing tomography surveys and 
CMB Planck data could be used to probe non-standard neutrino interactions and in particular the 
proposed new interaction with the CC. As it turns out the achievable explanatory power about the 
underlying interaction increases with the neutrino mass. 

In general, the combined projected sensitivity of lensing surveys and future CMB data is predicted to 
reach cr{Y^ rriy) > 0.03 eV for non-interacting neutrinos depending on the employed data sets and on 
the number of free parameters of the model (see Ref. [5] and references therein). According to solar 
and atmospheric neutrino oscillation experiments the neutrino mass squared differences are given by 
/^ml,^ = |Am^3| ~ 2.4 x IQ-^ eV^ and Am^^^^ = Amfa ^ 7.9 x 10"^ eV^ (e.g. Ref. [66]). Correspond- 
ingly, even if the lightest neutrino is massless, the most ambitious projects provide a 2a sensitivity to 
the minimal normal hierarchy with ^ ray = 0.06 eV and could thus distinguish it from the minimal 
inverted hierarchy corresponding to ^rrii, = 0.1 eV. Accordingly, non-standard neutrino interactions 
could be revealed by these future probes, if neutrino free streaming signatures are not observed (see 
e.g. [II])- Interpreted in terms of the considered scenario, this would indicate the sum of neutrino 
masses to be in the sub-eV range such that the corresponding non-standard relic neutrino density is 
reduced below the sensitivity of these experiments (cf. Sec. 15. 21 for an alternative, complementary way 
of testing the nature of the underlying interaction in the laboratory which is independent of the actual 
neutrino mass scale, since possible modifications to the neutrino distribution function are probed.) 

However, if X^m-i, = 1.2 -^ 2.7 eV as suggested by some members of the Heidelberg-Moscow collabora- 
tion [16], then it seems feasible for future cosmological probes to unambiguously reveal the proposed 
new neutrino interaction with the CC and to verify the Heidelberg-Moscow claim within this scenario 
as discussed in the following. 

Since CMB data and LSS surveys give us a snap shot of the universe at different times, they offer 
the possibility to compare the influence of neutrinos on cosmology before (or while) and after the 
interchange with the CC, respectively. 

If the sum of neutrino masses is in the super-eV range, the phase of interaction between the neutrinos 
and the CC overlaps with the epoch of recombination and thus during this time leads to a characteristic 
modification of the expansion rate as described above. Within a 7 parameter AMDM framework, CMB 
data is able to provide a bound on neutrino masses which is independent of LSS data due to the 
absence of a parameter degeneracy [5]. Considering that the projected sensitivity of Planck data alone 
is (t(^ rriu) < 0.48 eV for non-interacting neutrinos (see Ref. [5j for a review), it thus seems feasible to 
verify the Heidelberg-Moscow claim of X^mj^ > 1.2 eV for interacting neutrinos by CMB data alone. 

In addition, this result could be tested in a complementary way by future weak lensing or high redshift 
galaxy surveys which probe the relic neutrino background after the energy-momentum exchange with 
the CC has considerably reduced its energy density. This, however, implies uj = —1 = const, for the CC 
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Figure 4: The apparent neutrino mass as a function of redshift as given by Eq. ([2T]) for rrii^ = 2.2 
eV, nil, = 0.5 eV, rrii, = 5 x 10~^ eV and m,^ = 9 x 10~^ eV. Note that the ratio mapp/rrii, ~ 1/30 is 
independent of the neutrino mass for 1 + z < m/(9TQ). 



during structure formation and thus the absence of an apparent degeneracy between rrii^ and oj usually 
arising from a possible time-variation of w. Thus, it seems reasonable to compare to the projected 
sensitivity of future LSS probes within a 7 parameter AMDM framework. Since "^rrii, > 1.2 eV would 
correspond to Xl'^app > 0.04 eV, the most ambitious projects (see Refs. [Ml E]) promise to be able 
to verify the whole neutrino mass range indicated by part of the Heidelberg-Moscow collaboration for 
neutrinos interacting with the CC. In order to be proven right, the extent of the suppression of the 
power spectrum would require ^ wiapp > 0.04, but the corresponding scales of the power spectrum 
affected by neutrino free streaming would be > fcnr = 8 x lO^^/iMpc"^ according to Eq. ([27|l instead 
of > 2 X 10~^/iMpc~^ as expected from the standard relation in Eq. ([26]l . In summary, if the sum of 
neutrino masses is in the super-eV range, the prospects seem good to probe the considered scenario by 
future LSS surveys combined with future CMB data. In addition, in this case the absolute neutrino 
mass scale could be ascertained by identifying the predicted characteristic non-standard neutrino free 
streaming signatures on the power spectrum. 



5.2 Pauli-blocking efFects and interacting neutrinos 

Before we conclude this section, we would like to mention another interesting way of testing the consid- 
ered scenario in the laboratory which is complementary to cosmological measurements. Importantly, 
it turns out to be independent of the actual neutrino mass scale realized in nature. Namely, proposed 
experiments searching for Pauli-blocking effects [67l [68] resulting from the presence of the Cz^B are 
in principle sensitive to possible modifications of the neutrino distribution function / induced by the 
neutrino interaction with the CC as discussed in the following. For instance, in the experiment pro- 
posed in Ref. [68], the pair production of neutrinos with very low momenta p « is investigated, 
where the event rate is sensitive to the Pauli blocking factor (1 — /). In this context, a standard FD 
distribution /fd for the Cz^B would lead to a suppression factor (1 — /fd) — 1/2 in the event rate. In 
contrast, in our set-up, the distribution of the relic neutrinos is maximally degenerate with all states 
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occupied up to the Fermi momentum pp after the interaction with the CC has ceased, thus leading to a 
sharp-edged step function (1 — f{p)) = 6(p — pf) for the Pauh blocking term. Consequently, neutrinos 
cannot be produced with momenta below the current Fermi momentum pp = peTo/Te = yeTo, while 
the production of neutrinos with larger momenta is not suppressed at all. Compared to the FD case 
this implies a much sharper pair production threshold and in addition a shift in the threshold energy 
due to the non- vanishing Fermi momentum pp k, \.{) x 10~^ eV. 

6 Conclusions and Outlook 

The claimed evidence for neutrinoless double beta decay translates into a neutrino mass bound of 
^ m,y > 1.2 eV at 95% ^6] which is in tension with current neutrino mass bounds derived from cos- 
mology [6l [71 El m [TOl dU [III [TSj [m US]. In this work, we have shown that cosmological neutrino 
mass bounds can be relaxed and brought into agreement with the Heidelberg-Moscow claim, if a newly 
proposed neutrino interaction is taken into account which can be tested by future CMB and LSS 
measurements. It acts between relic neutrinos and the CC and arises from the zero-point energy con- 
tributions of the neutrino quantum fields. Since they induce a scale-dependence of the vacuum energy 
via renormalization group effects, the CC becomes time-dependent as long as the renormalization scale 
runs. In this case, an energy-momentum exchange between the relic neutrinos and the CC is implied 
through the Bianchi identity. 

We have studied in detail the time-evolution of the coupled system, in particular its consequences for 
the dynamics of the CC, the spectrum and the abundance of relic neutrinos and its impact on the 
interpretation of cosmological measurements. We have found that owing to the relative smallness of 
neutrino masses, the interaction becomes of dynamical influence at late times when neutrinos turn non- 
relativistic, Ty ~ rayj^. Moreover, it effectively ceases after a relatively short period, at a temperature 
also set by the neutrino mass, Ty ~ my/^. This decoupling behavior can be attributed to our particle 
specific choice for the renormalization scale. Since it was taken to be the average available neutrino 
energy which asymptotically runs to a value slightly above the neutrino mass, its evolution and thus 
the running of the CC effectively cease in the non-relativistic neutrino regime. By taking into account 
other fermions of higher masses, it was also shown that the corresponding vacuum contributions are 
becoming more and more important at late times, whereas they are completely sub-dominant in the 
early universe. During the cosmological evolution, the CC is therefore approaching a tracking regime 
at late times. 

Accordingly, in the relativistic neutrino regime after decoupling we have found the neutrino distribution 
to essentially maintain its standard Fermi-Dirac form, while the CC runs only logarithmically with 
time, but with a smaller value than today. However, when turning non-relativistic, the neutrinos were 
found to efficiently transfer energy-momentum to the CC, thus driving its energy density to its value 
as measured today. Since as a result, high neutrino momenta turned out to be removed from the 
neutrino spectrum while low momenta got enhanced, remarkably, the interaction was found to deform 
the neutrino momentum distribution into a maximally degenerate form. In comparison to Fermi-Dirac 
distributed neutrinos of equal mass, the mean neutrino momentum decreased by one order of magnitude 
implying a reduction of the neutrino energy density to 1/30 of its standard value after the interaction 
has ceased. Interpreted in terms of the standard relation between the neutrino energy density and mass, 
the non-relativistic neutrinos redshift as CDM, but appear lighter rriapp = 1/30 m,y. Accordingly, 
the presence of the non-standard neutrino interaction was found to considerably relax cosmological 
neutrino mass bounds. In particular, within the considered scenario current galaxy redshift surveys 
were shown not to be sensitive at all to neutrino masses below the upper bound ruy = 2.2 eV from 
tritium experiments ^U [20] and could thus be evaded. In addition, we have argued that the present 
comparatively mild neutrino mass bound from WMAP3 year data of ^ 2 ^ 3 eV at 2o" can be relaxed 
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in the presence of the proposed neutrino interaction with the CC. The reason primarily is that in 
comparison to standard ACDM cosmology, the characteristic integrated effects on the expansion rate 
turned out to be smaller than in the case of a AMDM model. 

However, we have proposed possible tests for the non-standard neutrino interactions with the CC 
by future cosmological probes as well as by laboratory based experiments. In the former case, the 
explanatory power on possible new physics turned out to increase with the neutrino mass. Since future 
weak lensing surveys probe a broad range of comoving wavenumbers at different redshifts, they provide 
unbiased, tomographic information about the matter power spectrum. Combined with Planck data the 
most ambitious projects promise a 2a sensitivity to the minimal value in case of a normal hierarchy 
(X] "^i' ~ 0-06 eV) [5] such that a non-observation of neutrino free streaming signatures would provide 
a hint for non-standard neutrino interactions. Interpreted in terms of the considered scenario this 
would mean that the interaction has reduced the neutrino energy density below the sensitivity of these 
future probes implying the sum of neutrino masses to be in the sub-eV range. 

In this case it might still be feasible for proposed laboratory based experiments to reveal the unique 
signatures of a possible interaction with the CC by directly probing the neutrino distribution function 
through Pauli-blo eking effects [671168] . 

If, on the other hand, the sum of neutrino masses is larger than ^ my > 1.2 eV as suggested by part 
of the Heidelberg-Moscow collaboration, we found future LSS probes also to be able to unambiguously 
reveal the characteristic free streaming signatures of neutrinos having interacted with the CC. Namely, 
Fermi-Dirac distributed neutrinos can free stream larger distances the lighter they are. However, as 
a unique signature of the energy-exchange with the CC, neutrinos become slower than in the absence 
of the interaction and thus free stream a shorter distance, while at the same time they appear lighter 
due to their reduced abundance. 

In addition, it seems likely that future Planck data alone could provide a complementary test for 
the Heidelberg-Moscow result within the considered scenario, since in this case the interacting phase 
overlaps with the recombination epoch. Considering that the projected sensitivity in the absence of 
the neutrino-interaction is a{'^my) < 0.48 eV (see e.g. Ref. [5| and references therein), it thus seems 
feasible by Planck data to trace the characteristic modifications to the expansion rate resulting from 
the neutrino - CC interaction. 

Let us in the following give an outlook on possible extensions of the considered scenario. Since theories 
beyond the standard model predict other relics as warm or cold dark matter candidates, e.g. axions, 
gravitinos or supersymmetric partners of standard model fermions, their effects on the running of 
the CC relevant at energy scales of the order of their masses, could consistently be studied within the 
proposed renormalization scheme. Such an extension of the considered scenario is especially interesting, 
because the masses of bosonic species enter the RGE for the CC with an opposite sign as fermions. 
This might indicate that while the interaction with the CC generically causes the fermionic species to 
lose energy-momentum, in turn, the bosonic species would gain energy-momentum. Thus, presumably 
this would tend to increase the momenta of the bosonic distribution, making the bosonic dark matter 
candidate hotter, while the energy density in the CC would decrease. However, in order to be consistent 
with observations, within this scenario strong restriction on the particle masses are expected to arise 
or their existence could turn out to be excluded after all. 

After the posting of the first version of this paper, Ref. [69] appeared. The corresponding authors 
claim that in the work on the RGE running of the CC [39i i40j part of the one-loop corrections to the 
CC have been missed out, namely the explicit // dependent part of the full one-loop effective action. 
It is furthermore argued that on including the missing part the CC does not formally run with the 
renormalisation scale /i. In our opinion this claim is misleading, since it is precisely this allegedly 
missed out one-loop correction which gives rise to the non-zero beta-function for the CC in Eq. ([T]) 
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(determined by the partial /i— derivative). Moreover, in the framework of an RG improvement within 
the effective field theory approach taken in our work, it is as usual the (non-zero) beta-function which 
describes the RG evolution. Therefore, we do not see the appropriateness of our RG analysis of the 
CC invalidated by the above paper. 
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